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ON  SUBSONIC  FLOW  OF  A COMPRESSIBLE  FLUID 


P.  R.  Garabedian  and  M.  Schiffer 


The  results  presented  here  were  motivated  by  a desire  to  give  a simple 
treatment  of  the  non-linear  elliptic  partial  differential  equation  governing 
the  steady  irrotational  subsonic  flow  of  an  ideal  compressible  fluid.  For 
two  independent  space  variables  x and  y,  the  stream  function  4*  of  such  a 
flow  satisfies  an  equation 

(1)  4’y)’4'yi"  o 

where  is  an  analytic,  increasing,  convex  function  cf 

2 2 l/?.  ^ 

q-  (4  + ^ ) whose  explicit  form  depends  on  the  equation  of  3tate  of  the 
x y 

fluid  in  question.  Our  analysis  of  (l)  will  be  based  on  the  fact  that  it  is 
the  Euler-Lagrange  equation  for  the  double  integral  variational  problem 


fj'  y^)dxdy  " minimum. 


We  shall  introduce  several  devices  for  analyzing  (2)  which  prove  to  be 
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particularly  successful  for  the  case  F-  (l+4  + 4 ) of  the  Plateau  problem. 

x y 

Shiffraan  [5]  has  given  a proof  of  the  existence  of  subsonic  compressible 
flows  based  on  (2).  his  work  is  part  of  an  extensive  literature  cn  the  calculus 
of  variations  and  on  non-linear  elliptic  partial  differential  equations  of  which 
the  contributions  by  Haar  [2]  and  Rado  [4]  come  closest  to  the  point  of  view 
of  our  paper.  The  deduction  of  a priori  estimates  on  the  derivatives  of  the 
solution  4 of  (2)  and  the  discussion  of  the  analyticity  of  4 are  key  develop- 
ments in  the  theory,  and  it  is  in  these  two  directions  that  our  analysis  applies. 
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In  the  second  section  of  the  paper,  we  study  the  minimum  problem  (2)  by 
the  method  of  interior  variation  [1].  This  leads  in  a natural  way  to  an  integral 
H for  which  we  can  derive  a second  order  partial  differential  equation  from  the 
existence  of  th9  first  derivatives  of  In  the  case  of  the  Plateau  problem, 

H satisfies  an  elliptic  Monge- Ampere  equation  which  ties  in  with  Rad6's  proof 
of  the  analyticity  of  a minimal  surface. 

As  a preparation  for  the  application  of  interior  variations,  we  take  up 
in  the  third  section  a construction  based  on  symmetrization  which  yields  for 
the  Plateau  problem  a minimal  sequence  satisfying  a uniform  Lipschitz  condition. 
For  this  construction,  our  assumption  on  the  boundary  data  is  more  general  than 
that  usually  required  for  the  analogous  conclusion  using  the  three-point 
condition,  and  we  are  therefore  able  to  discuss  the  Plateau  problem  in  non- 
parametric  form  for  a domain  which  need  not  be  convex.  Furthermore,  the 
symmetrization  method  applies  equally  well  for  any  finite  number  of  indepen- 
dent variables . 

2.  iaterlQE ..YarlaUQh.aufl  ..ths..  late  seal  H. 

In  a plane  domain  D with  boundary  curves  C.  let  'f  be  a solution  of  the 

minimum  problem  (2),  with,  for  example,  prescribed  boundary  values  on  the 

curves  C.  We  shall  assume  merely  that  41  satisfies  a Lipschitz  condition  in  D, 

so  that  the  first  derivatives  41  and  4*  exist  almost  everywhere  and  are 
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bounded.  Instead  of  trying  to  derive  the  Euler- La grange  equation  (l)  by  the 
classical  method  of  varying  SP,  we  study  (2)  in  this  section  by  performing 
infinitesimal  transformations  on  the  independent  variables  x and  y and  by 
considering  the  shift  of  ^ thus  generated. 

We  let  f be  a continuous  complex- valued  function  in  D,  with  piece-wise 
continuous  first  derivatives,  which  vanishes  on  C.  It  is  convenient  to  use  the 
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complex  notation  2f  - f -if  . 2f_- -Mf..,  z • x+iy,  for  derivatives.  For 

Z X „ z A ¥ 

snail  values  of  the  complex  parameter  € , the  transformation 


z - z ♦ € f 


-*  * + 4 * 

z ■ x ♦ iy 


performs  a one-to-one  mapping  of  D onto  itself.  We  define  a new  function  4^ 
in  D by  th9  formula 

U)  y (x*,y*)  - 4^(x,y)  , 

and  we  compare  the  values  of  the  integral  in  (2)  associated  with  the  two 

* 

neighboring  functions  4*  and  ^ . 

Clearly 

(5)  Jj  »ur:  KJix-  4/ 

whence  by  elementary  calculation 


p(44/*4i*)dx  <*y  “ j'j" r(4yz  4L)dxdy 


-2  Re  [ |C(F-  q2F')f  -4F'^  f.ldxdy  i + 0(|e|2) 


where  q » 4 4^  4L  and  where  F and  F'  stand  for  F(q2)  and  F'(q  ).  In  the 

z z 

usual  fashion,  we  conclude  from  (6),  from  the  extremal  property  (2)  of 

and  from  the  arbitrary  nature  of  € that 

(7)  JJ[(F-  qV)fz- 4^'H^  ^Idxdy  - 0 

for  every  continuous  piece-wise  continuously  differentiable  function  f in  D 
which  vanishes  on  C. 

Let O be  a closed  subdomain  of  2,  let  co  be  a continuously  differentiable 
function  in  D which  is  identically  1 in  Cl  and  which  vanishes  on  C,  let  t be  an 
interior  point  of  Cl,  and  let  p>0  be  so  sma]]  that  the  circle  |z-t|^  ^ is 
contained  in  the  interior  of  O . We  define 


- 4 - 


! . 


! • 

: 

l 


i 

t 

t 

t 


(3) 


f . 

e2 


for  Iz-tl  p and  we  define 

60 


(9) 


f - 


z-  t 


in  D for  |z-t|  We  can  substitute  this  function  f into  (?)  to  obtain 

(10)  -Aj  IT  F'^dx  dy  - - JX  — dxdy  ♦ A(t)  , 

P2  ^ 2 Cl-K  (z-  t)2 

where  K denotes  the  circle  [z-t|^r  p and  where  A(t)  is  the  analytic  function 

(11) 


A(t)  - (P-  qV)[^  - — ^)dxdy 

(z-  t) 


D-0 


in  O . 


Letting  p->0,  we  find  almost  everywhere  in  Ci 

(12)  - - jj 


n 


gr~a— -■  dxdy  + A(t) 
(z-t)* 


where  the  integral  on  the  right  is  to  be  interpreted  in  the  sense  of  the  Cauchy 

2 

principal  value.  In  the  case  of  the  Dirichlet  problem,  F =.q  and  the  formula 

2 

(12)  shows  immediately  that  4^  is  an  analytic  function.  Thus  we  obtain  in 
the  simplest  situation,  corresponding  to  an  incompressible  fluid,  a quite 
elegant  proof  that  the  solution  of  the  minimum  problem  (2)  is  a regular  function. 

In  order  to  study  the  general  non-linear  problem,  we  introduce  the  integral 

(13) 


H--^r^J(F-  q2F'  ) log  [ z-t  | dxdy  + B , 


where  B i3  a real  harmonic  function  in  O.  3uch  that  TTB.^  - -A(  t) . By  (12) 
and  by  standard  lemmas  on  the  second  derivatives  of  a logarithmic  potential, 
we  find  almost  everywhere  in  O. 

(H>  H 


-F-qV  , 


(15) 


tl 
®tt  - 


H - -4P'  vl/2 

tt  t 


E 
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This  give 3 in  turn 

(16)  H. +H  - H2  *2q2FF'-F2 

tt  tl 

ffe  can  eliminate  q from  the  equations  (14)  and  (16)  to  obtain  for  the  real 
function  H a partial  differential  equation  of  the  form 

(17)  , 

where  Q is  a real  analytic  function  of  its  real  argument  which  is  completely 
determined  by  F.  For  the  most  general  function  F corresponding  to  an  arbitrary 
equation  of  state,  (17)  is  a non-linear  equation  equivalent  to  (l)  which 
involves  only  very  special  combinations  of  second  derivatives  of  the  auxiliary 
function  H.  The  significance  of  this  second  order  partial  differential  equation 
is  that  its  derivation  requires  only  a Lipschitz  condition  on  the  stream 
function  MA 

The  form  of  (17)  suggests  finding  those  integrands  F for  which  it  reduces 
to  the  Monge-Amp&re  equation 

AS) 

This  reduction  takes  place,  according  to  (14)  and  (l6),  when  F 3ati3fie3  the 
ordinary  differential  equation 

(19)  2q2FF' - F2  - 2A(qV- F) 

for  a suitable  value  of  th9  constant  A . We  check  immediately  that  ( 19 ) has 
the  general  solution  

(20)  F - A ♦ ^ A2*  |J?q2' 

whence  the  Ifcnge-Amp^re  equation  ( 18 ) is  3een  to  correspond  to  the  case  in 
which  (2)  is  the  Plateau  problem. 

We  arrive  in  this  way  at  a proof  of  the  analyticity  of  a minimal  surface. 
For  we  can  apply  the  Legendre  transformation 
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(2l)  h+H  - ux+vy  , u -*  , v-Hy 

to  (18)  in  order  to  obtain  the  linear  elliptic  quation 


(22) 


h ♦ h - 1 
uu  w 


i'he  solution  h of  (22)  is  evidently  an  analytic  function  of  u and  v,  whence  H 
is  an  analytic  function  of  x and  y,  and  by  (14)  and  (15)  we  can  conclude  that 
the  function  is  also  analytic  in  x and  y.  The  Poisson  equation  (22)  yields, 
furthermore,  a procedure  for  constructing  flow  patterns  explicitly  according  to 
the  formulas  of  this  section  in  the  case  where  the  equation  of  state  of  the 
fluid  leads  to  an  integrand  of  the  type  (20) . 

3.  Symmetrlzatlon  and  the  llpschitz  condition. 


The  analysis  of  the  preceding  section  exploited  formal  manipulations  in 
order  to  demonstrate,  in  certain  important  special  cases,  the  analyticity  of 
solutions  of  (2).  In  this  section,  we  complete  our  discussion  of  the  Plateau 
problem  by  constructing  a minimal  sequence  which  fulfills  a uniform  Lipschitz 
condition. 

Along  the  curves  C bounding  the  domain  D,  we  assign  boundary  values  which 
generate  in  space  a system  of  smooth  closed  curves  P.  The  problem  is  to  span 
through  the  curves  P a non- parametric  minimal  surface  over  the  domain  D.  For 
this  problem  of  the  type  (2),  we  can  clearly  find  a minimal  sequence  of  surfaces 
whose  areas  approach  the  minimum  value  in  question.  There  is  no  loss  of 
generality  if  we  assume  that  each  of  these  surfaces  lies  in  the  convex  hull 
of  the  system  of  curves  P,  since  when  this  is  not  true  for  a given  surface, 
we  can  diminish  the  area  of  the  surface  by  replacing  portions  of  it  by  sections 
of  planes  tangent  to  the  convex  hull  of  P.  Furthermore,  it  is  permissible  to 
diminish  the  area  of  any  of  the  surfaces  of  the  minimal  sequence  in  a similar 
manner  by  replacing  portions  of  them  by  simply- connected  sections  of  catenoids, 


i 

t 
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or  of  other  specific  known  minimal  surfaces,  which  do  not  intersect  P. 

Thus  we  may  suppose  that  all  members  of  our  minimal  sequence  lie  in  the 
largest  closed  region  E projecting  onto  D which  contains  the  curves  P,  but 
which  does  not  intersect  in  a simply- connected  surface  element  any  plane  or 
catenoid  not  meeting  P.  This  latter  condition  means  that  E cannot  be 
diminished  infinitesimally  by  cutting  off  a volume  element  with  a plane  or 
catenoid  which  does  not  intersect  P. 

We  now  make  the  assumption  on  the  boundary  curves  P that  the  closed 
set  E is  ao  situated  that,  for  some  G>  0,  all  lines  making  an  angle  smaller 
than  0 with  the  normal  to  the  plane  domain  D and  intersecting  P have  only 
one  point  in  common  with  E.  This  assumption  restricts  the  curvature  of  P, 
but  does  not  imply  that  D is  a convex  domain,  3ince  in  some  cases  D can  even 
be  multiply- connected. 

We  consider  any  rectangular  coordinate  system  in  which  the  z-axis  makes 
an  angle  les3  than  6 with  the  normal  to  the  plane  of  D.  In  this  coordinate 
system,  any  element  of  our  minimal  sequence  has  a non- parametric  representation 
z»z(x,y)  which  may  be  multiple-valued,  with  branches  z = z^(x,y),z  - Z£(x,y) , . . . , 
z ■ z^+i^y)  > z^£si  z2m-fl  * sy^^rize  such  a surface  by 

replacing  it  with  the  surface  whose  non- parametric  representation  is  the 
s ingle- valued  expression 

2m+l 

(23)  z - y v (-l)k*1  z.(x,y) 

k- 1 K 

It  follows  from  the  basic  results  of  Steiner,  or,  more  directly,  from 
Minkowski's  inequality  [3],  that  the  syrametrization  process  (23)  diminishes 
or  leaves  unchanged  the  area  of  the  surface.  Furthermore,  if  symmetrization 
in  one  such  coordinate  system  is  followed  by  symmetrization  in  another,  the 
resulting  surface  still  has  a single-valued  non- parametric  representation  in 
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the  first  coordinate  system.  This  can  be  checked  by  making  en  affine  transfor- 
mation such  that  the  directions  of  the  two  synmetrizations  become  perpendicular, 
a case  in  which  the  result  is  evident.  Finally,  and  nost  important  of  all,  the 
symnetri zation  procedure  (23)  does  not  alter  the  boundary  curves  P of  any  of 
our  surfaces,  since  the  surfaces  lie  in  the  set  E,  which  has  the  property  that 
if  any  line  parallel  to  one  of  our  z-axes  intersects  P,  then  it  intersects  P 
and  E in  only  one  single  point. 

These  considerations  show  that  we  can  assume  without  loss  of  generality 
that  each  surface  of  our  minima]  sequence  has  a single- valued  non- parametric 
representation  in  every  coordinate  system  whose  z-axis  is  inclined  at  an  angle 
less  than  0 with  the  normal  to  the  plane  of  D.  But  if  we  denote  by  z » ^(xjy) 
the  representation  of  such  a surface  in  a coordinate  system  such  that  D lies 
in  the  (x,y)-plane,  then  this  result  implies  that  "41  satisfies  the  Lipschitz 
condition 

(24)  |'V(x2,y2)- m[(x2-x1)2+  (y2-y1)2l  ^ , 

with  M-  cot  0.  Thus  we  obtain  a minimal  sequence  of  surfaces  satisfying 

the  Lipschitz  condition  (24),  end  from  the  lewer  semi- continuity  of  the 

area  integral  [2]  we  can  deduce  the  existence  of  a solution  of  (2)  satisfying 

2 l/2 

the  same  Lipschitz  condition,  for  F-  (l  + q ) 

A combination  of  the  techniques  of  this  section  and  of  the  previous 
section  yields  a solution  of  the  Plateau  problem  in  nor.-peranetric  form  for 
a domain  D which  is  not  necessarily  convex,  but  for  a system  of  smooth  boundary 
curves  P satisfying  the  above  geometrical  condition  relative  U/  certain 
planes,  catenoids , and  projections.  Oar  method  for  developing  such  a condition 
is  a generalization  to  non-linear  elliptic  equations  of  the  majorization 
principles  of  the  theory  of  linear  elliptic  partial  differential  equations 
based  on  the  maximum  principle. 
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A point  of  interest  in  the  symmetrization  construction  presented  in 
this  section  is  that  it  yields  estimates  of  precisely  the  same  nature  in  space 
of  any  number  of  dimensions.  In  particular,  we  can  treat  the  variational 
problem 

(S5)  ‘ dxdydz  * miniicuE, 

which  has  applications  to  the  study  of  three-dimensional  flow  of  a 
Karman-Tsien  gas. 
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